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1 Introduction 

The question of the composition of the ultimate and densest phase of matter real¬ 
ized in nature is a fundamental problem in physics. Approaches to a partial answer 
have arisen from different directions. From the theoretical point of view, in the early 
universe the condensation from a quark phase to a nucleon phase has already been 
studied using cosmological arguments [1]. On the other side, at RHIC experimental 
efforts from the heavy-ion community trying to produce the so-called quark-gluon 
plasma have yield even more questions on the nature of this phase of matter [2] . The 
possibility of quark deconhnement seems possible in a very dense and/or high temper¬ 
ature environment according to the quantum chromo dynamics (QCD) temperature- 
baryonic density phase boundary T{pb)- In this context, it is especially interesting 
that the conditions believed to exist inside dense cold stars are those matching the 
low-temperature, high-density region of the before-mentioned phase space. Neutron 
stars (NS) may constitute, in this way, possible sites where nucleon and quark phases 
can coexist and serve as a test-bench at the extreme for fundamental interactions. 

Besides standard-model baryonic matter, the presence of another type of matter 
present in the universe, dark matter (DM), may add interesting and new insights to 
current standard scenarios. This type of matter makes more than 80% of the current 
matter content of our universe and has been conhrmed by precision measurements 
from a range of areas, for example, rotation curves in our galaxy, such as those recently 
presented in [3]. The intimate nature of the dark matter remains to be determined, 
i.e. its fermionic or bosonic nature or its couplings with the rest of the standard- 
model particles. There is a huge world-wide effort to constrain its properties either in 
a direct, indirect or collider search. These missing pieces of information will be key 
in revealing the possible observable signatures of its existence. 

In this contribution we elaborate on the consequences of the presence of DM in 
a dense nuclear matter star environment focusing on indirect effects, namely that 
of bubble nucleation. Let us mention that this same output is already considered 
as a direct search product in experimental settings with detected events consisting of 
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Figure 1: 2L PICO bubble chamber (left panel) and MOSCAB experiment based on 
the geyser technique (right panel). Figures courtesy of both collaborations. 


recoiling of a nuclear target after a weakly interacting particle hits on it. Experiments 
along these lines include PICO [1], SIMPLE [5j or MOSCAB [6], see Fig.([^. 

The consideration of dense star settings requires the scaling of the same type of 
physics used in bubble chambers or superheated droplet detectors in extreme envi¬ 
ronments with differences of several orders of magnitude in density and temperature. 
This is a real theoretical challenge and requires state-of-the-art microscopic calcula¬ 
tions at the quantum level. Although quantum bubble nucleation has already been 
discussed theoretically in the literature [7] in the context of quantum fluctuations of 
temperature or seeding |S], it remains to be tested on earth from the experimental 
point of view. The key to this process is the energy deposition that could lead con- 
hned quarks out of the nucleon bag, overcoming a ~ (0.1 — 1) GeV potential energy 
barrier. As we will argue the presence of external energy sparks due to DM particle 
self-annihilations could indeed drive the deconhnement transition. 


2 Dark matter as a nucleation trigger source 

In the classical bubble nucleation theory of Seitz [9] an initial localized energy de¬ 
position over a threshold value, Eg, in a dense medium, is the mechanism to drive 
microscopic bubble formation. Further modellization allows that a typical ionization 
reaction involving complex processes of atomic collision-cascades is replaced by an 
abrupt temperature rise in an inhnitesimal cylindrical volume around the ion trajec- 
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Figure 2: Content of the current universe. Dark matter amounts to 23%. Figure from 
Chandra observatory website at http://chandra.harvard.edu 


tory at the time-of-passage. Activated particles can populate a metastable state with 
energy E with a probability ~ Q-E/kgT^ j£ high temperature endures only for a 
short time, the higher states have only a small chance to return to the ground state, 
due to the potential-energy threshold associated. 

In a NS core, the temperature and density conditions are those of a degenerate 
Fermi system. Typical nucleon energies are those of the Fermi level, fulhlling the 
relationship E/ksT ~ Ep/ksT S> 1. For such a system mass density is thought to 
be several times that of nuclear saturation, po — 2.8 x lO^^g/cm^, and it is highly 
improbable to excite fluctuations of 5T ~ 10 MeV estimated to be necessary for the 
quark deconhnement [TO] . 

In this pictured scenario an external agent is needed to perturbe the nucleon 
degeneracy in order to allow the system to free some of the quark content. At this 
point we must include the dark matter fraction. It amounts to 23% of the current 
content in our universe, 8 times larger than the standard-model particle content, see a 
cartoon in Fig. (|^. Let us consider a Majorana DM component, electrically neutral, 
capable of self-annihilating and massive enough to be gravitationally accreted from an 
existing galactic DM distribution onto a NS. At a given galactocentric location close 
to the sun the local DM density ~ 0.3GeV/cm^ is taken as reference. Particle 
candidates of this kind are indeed in the literature na, one of the most popular, the 
neutralino, belongs to the supersymmetric class. For a mass range ~ (10 — 10"^) 
GeV/c^ and inside the NS, we will take a (spin independent) scattering cross section 
with nucleons around the sensitivity values of current direct detection experiments. 
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a^jsf ~ cm^. This choice allows building a stabilized thermalized fraction of DM 
aside the standard baryonic component over time. Using this fact one can write the 
DM particle number density [13] as 

= no,^e~(^) , ( 1 ) 

where rth is the thermal radius rth = \J 2 ^ 0 central DM density 
value normalized to the global DM population inside the star of mass M and radius 
at a given time, N^(t). 

We take DM self-annihilation as the only process removing DM from the star. 
The X"Population number is obtained from the solution of an ordinary differential 
equation — CaN^- It explicitly includes competing processes from capture 

and annihilation with rates and Ca respectively [T3| . The capture rate is given by 

^ 3 X 10^^ ( M \ f R \ / lTeV \ ( 

^ /(f) Vb5M0y VlOkmy \ J 1^0.3^ 

with / (f) = 1 - 0.4 (^) a redshift correction factor. Following [I5] we 

do not consider a reduction of the cross-section since we use values well above the 
geometrical cross-section. It is important to note that this dependence sets a limit 
on the intrinsic capability of NS to accumulate a critical amount of DM and possibly 
serve as a test-bench for DM properties. For the annihilation rate we write 



a 



(r, T){aav) dV. 


Then we can hnally write 


NPt) 



+ coth-‘ 

T 



(3) 

(4) 


with = ^jC^Ca the inverse relaxation time to achieve equilibrium and iV^(tcoi) 
is the number of DM particles inside the progenitor core at the time of the collapse 
(NS birth). This population is essentially inherited from the progenitor star in its 
lifetime. 


3 Induced bubble nucleation 

The possibility to nucleate a quantum bubble of quark matter will be driven by the 
ability of the DM self-annihilation process to effectively inject the spark of energy 
needed to activate quarks conhned in the nucleon bag. Note that the ultimate goal of 
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this process could be the formation of a bubble, in a similar way to thermal boiling. 
Since there are however some slight differences among both processes, we will however 
assume here that the energy required is the same for the sake of brevity. 

The so-called free energy of nucleation is approximately H! given by F ~ 

Pq) + StijR where the two terms account for the pressure and curvature energy con¬ 
tributions, respectively. In order to produce spherical bubbles the minimum of this 
quantity shows a critical radius, namely = 0 implies 


where AP = Pq — P^ and Pq (P„) is the quark (nucleon) pressure and stability is 
granted as the contained baryonic number A ~ ^nR^nn > ^min and Amin ^ 10 [TB] . 
We will consider here, for simplicity, that the system is pure neutron matter. 

In the approximation of very light quarks (i.e u,d quarks), fUq ~ 0 and creation 
of neutrally charged bubbles this critical work is given by 


F{Rc) = IV, 


IbTT / 27^ 

irV AP' 


( 6 ) 


For a two-flavor ud-quark system in the MIT bag model (and B > (145 MeV)"^) [7]) is 

4 { 1 ‘2\^ /I 

given hy Pq= 'Yh ~ ^ and assuming a pure n system P^ ~ rnost of 

i=u,d 

_ _ 2 

the pressure will effectively be provided by neutrons in the n fluid [7]. 7 = Y is 

i=u,d 

the curvature coefficient and pj (p^) is the quark (nucleon) chemical potential related 
to the Fermi momentum of the degenerate system pi = ppi {pn = + p%^). 

Electrical charge neutrality requires for the ud matter number densities na = and 
nn = {nu + Urf)/3 with Ui = ^ in the light quark massless limit. 


3.1 Activation of quarks 

In this subsection we grossly estimate the injection of energy from hnal state photonic 
channels coming from the DM self-annihilation. We will discard the many-body 
blocking in nuclear phase space, estimated to be a minor effect at this range of 
energies. Additional channels should be considered carefully, although large masses 
> 1 GeV/c^ allow a safe semiclassical treatment of the DM component. The rate 
of external energy (photon channels) from the annihilations can be written as 

dE f 

— = j EQ{E,r,T)dEdV, (7) 
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where Q{E,r) is a function defined as a sum for each possible annihilation channel 
as 

Q{E, r,T) = Y^ T) ~ nj(r, T) (§) 

i 

dN'^ 

is the photon spectrum for the i-th annihilation channel. Let us note that the 
low-energy region of the spectrum will not be capable of activating the quarks if 
E < Eq, then, an efficiency of the process, /^ < 1, is expected. However, in terms of 
mass fractions for example, even a /« = 0.1% for a ~ TeV particle would allow 
the triggering. It is worth mentioning at this point that, in principle, also decay 
modes could serve as an activation method [H] , however we will restrict ourselves to 
annihilation in this work. 



r/R 



r/R 


Figure 3: (Left panel) Internal temperature radial profile in a NS with a self- anni¬ 
hilating DM component with = 10 TeV/c^ at a t = 0.1 yr since collapse. (Right 
panel) Inner core range < 1 km [20] • See text for details. 

Fig. ([^ panel on the left shows internal temperature vs. radius in the star. The 
local heating is compared in a standard (MURCA neutrinos [12] and surface photon 
emission) scenario (solid line) and that due to DM annihilation for a = 10 TeV/c^ 
(dashed line). On the right panel the local temperature radial profile is zoomed to 
the inner ~1 km. Inner early heating could produce effective boiling. Details are 
given in |20|. An arbitrary flat profile with T = 10® K is supposed at initial time 
and a value {(Tav) ~ 3 x 10“®® cm^/s for the velocity averaged DM annihilation cross 
section is taken. 

The time averaged energy density coming from this DM process can be estimated 
as 

1 rto+At p 

{udm)m -Xtj, 7 ( 9 ) 
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This is to be compared to the energies quoted to produce bubbles DU ; ^bub — 
Wc/Vd ~ 5.4 X lO^^erg/cm^, or in more familiar nuclear units Ubub ~ 0.1 GeV/fm^. 
The spatial spread is Ar ~ ^ ^ ^ ^2 — 2 fm, ensuring the local char¬ 

acter of the energy release. In heavy-ion collisions the estimated energy density for 
deconhnement is m > 5 GeV/fm^ at low pressures [21] or hve times smaller in lattice 
QGD calculations izg. It seems that a high pressure environment could catalyze the 
reaction, although further refinement of the models is still needed. 


4 Bubble associated instabilities 

In this scenario the number of nucleated stable bubbles created can be accumulated 
in the central region over a period of time. To produce an observable effect we expect 
that after a critical number of bubbles is created the coalescence of bubbles in the 
dense nucleon system may induce further macroscopic structure changes in the star. 
The quoted number for this to happen is however not clear at this point. Some 
authors like Harko and collaborators m state that one single bubble Nq = 1 will be 
responsible for driving a NS to quark star (QS) transition. This is due to a bubble 
dynamical mechanism allowing it to grow to become comparable to the star size. 
More in general, the relationship to analyze is given by 

,,, f dNbuhdE 

" J 

In this contribution, we will be conservative and will study the consequences of driving 
a moderate mechanical instability in the central core of the star. If this was the case 
the equation of state of baryonic matter P(T, p) would allow pressure oscillations that 
could disrupt the star producing a conversion into a more compact object. Whether 
the conversion is partial or complete is not yet completely understood. 

As a consequence of the triggering, the creation of a number of bubbles ^Wub 
could cause a perturbation of the hydrostatic pressure value 

^Abub, (11) 

0 

with respect to the standard value. The derivative coefficient is related to regular 
nuclear matter compressibility (in the absence of bubbles). To create an instability, 
this change should be such that the star can not adjust adiabatically to it so that a 
catastrophic event is driven. In order to quantify the effect, let us consider an old 
and cold NS with T ~ 10^ K and pn = 5po for a DM particle with ~ 1 TeV. In 
principle, a standard variation in the number of nucleons depleted from the thermal 
volume nucleon sea, ■\/Ath — 10 ^^, could provide a large enough variation in pressure 
to trigger transitioning. Then the number of bubbles to create this perturbation is 


SP ~ 


dP 


dN, 


bub 
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No ~ VA.hM miTi — 10^°. From eq.([^ we can see than there is a comparable captnre 
rate of DM particles into the star. At late times t » t eqnilibrinm can be achieved 
and an efficient mechanism conld indeed create snch a nnmber of babbles, precipitat¬ 
ing the macroscopic conversion of the nnclear star into a (hybrid) QS. The conver¬ 
sion itself has been extensively stndied in the literatnre see for example [m ESI ED 
where emission of gamma ray bnrst is likely to follow the rapid (ms timescale) event. 
Since agile methods of gamma-ray detection are increasingly available, characteristic 
featnres in the photon emission predicted in these events (sharp signal) conld help 
disentangle the formation of these compact objects and go deeper into the possible 
formation mechanisms in the fntnre. 
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